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Abstract: A persistent and fascinating problem at the high energy colliders are jets. 
Often trying to observe physics underlying the hard interactions at colliders requires ex¬ 
perimental cuts in phase space, defining several jet or beam regions. QCD being a gauge 
theory that readily decays into infra-red modes, correlations between jet regions is almost 
inevitable, spoiling the predictivity of fixed order QCD calculations. One is faced with 
the task of calculating the evolution of a reduced density matrix, where successively less 
energetic (jet) regions are integrated out, to gain control of the calculation. I relate the 
decay rates governing the flow into the IR to an effective field theory expansion in soft jets, 
allowing a systematic and resummed calculation of these rates, while further relating them 
to physically observable features of the QCD cascade. To demonstrate the utility of the 
soft jet expansion, I present a factorization theorem for a soft subjet collinearly splitting in 
and out of a parent fat jet. Using the resummation properties of this factorization theorem, 
I elucidate the structure of the subleading non-global logs (encoding the jet correlations) 
in the hemisphere jet mass distribution, as well as give a collinear improvement of the 
leading order resummation equation, the BMS equation. I compare to other approaches 
to subleading resummation of NGLs, and find the collinear improvement of the leading 
order equation removes the need for kinematic-dependent corrections in the IR averaging 
procedure of the reduced density matrix, so that no further large logs can be generated in 
the IR. Finally I end with speculation about connections with collinear improvements of 
the NLO B-JIMWLK hierarchy for small-x resummation. 
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1 Introduction 

Non-global logarithms (NGLs) pose an fascinating problem for resummation [1], and can 
have a large phenomenological impact [2, 3]. These are logarithms associated with soft 
radiation which splits and lands into different angular phase space regions, with an energy 
hierarchy between the soft daughters. Physically one pictures as an initial condition mul¬ 
tiple hard jets at wide fixed angles in the event, whose masses are all much smaller than 
the initial hard collision. However, the masses themselves exhibit a hierarchy of scales. 
The resummation is difficult, since the soft radiation is sensitive to the pattern of previous 
soft emissions. This is in contrast to global logarithms, the logarithms of the jet mass to 
the hard scale, which are associated with virtual corrections of QCD amplitudes. For the 
global corrections, no parton becomes a physical state between the hard interaction scale, 
and the jet mass scale, so that while the cancellation of IR divergences between real and 
virtual corrections still occurs, a large log is left over. The exponentiation of these global 
logarithms is now a standard technique, whether using factorization theorems [5-13], or 
QCD coherence arguments [14-18], to the point that it can be essentially automated [19- 
22] in both approaches. Moreover, the exponentiation is truly an explicit exponentiation, 
at least in an appropriate conjugate space for additive observables. However, the fixed or¬ 
der pattern of NGLs as calculated in [23] and [24] do not exhibit any such straightforward 
pattern. 

Dealing with the resummation of NGLs, four approaches have appeared in the lit¬ 
erature: avoiding them entirely with clever observables [25-27], a large-A^c Monte Carlo 
technique [1], non-linear evolution equations^ [28-35]^, and an expansion in soft sub-jets 
[38-42]. This last approach can be grounded in the structure of factorization theorems, as 
derived in the framework of soft collinear effective field theory (SCET) [43-45], promising 
a straightforward way to calculate to higher orders in Ug for the anomalous dimensions in 
the resummation. Perhaps even more importantly. Ref. [42] organized the phase space of 
emissions into parametrically separated regions where distinct (now global) resummations 
apply. This is accomplished by using a complete set of infra-red and collinear safe observ¬ 
ables, for example see Refs. [46, 47], to define multi-differential cross-sections with sufficient 
number of observables to distinguish the soft and collinear limits of multiple pronged jet 
structures [48-52]. Resumming these cross-sections, and expressing the more inclusive 

^For NGLs, this evolution equation is called the BMS equation after its inventors, Banfi, Marchesini, and 
Smye. However, it fits into a universality class of quantum master equations, see Sec. 11 for a discussion 
on the universality of IR evolution equations. 

^For progress towards an EFT interpretation of these evolution equations, see Ref. [36], which incor¬ 
porates an object akin to the trace of the post-evolution reduced density matrix of Ref. [37] into their jet 
factorization theorem. 
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cross-sections as a marginalization over the more exclusive ones, the perturbation series is 
re-organized as summing over jets (resummed partons), rather than fixed order partons. 
Interesting emergent behavior, like the buffer region [53], manifests itself for these dressed 
emissions. Since the observables are complete up-to a resolution scale (which terminates 
the number of jets included in the calculation), one can also check whether all phase space 
regions are accounted for by examining all possible relative scalings of the observables. 
This has important implications for the subleading NGLs, since the fixed order NGLs at 
in QGD are sensitive to a distinct phase space region, where a soft subjet is collinearly 
splitting along the jet boundary. This was already noted in Ref. [42] . Here I shall worked 
out the soft jet factorization in this region, and its implications for resumming NGLs. 

The organization of this paper is as follows: I set out the phase space region for a soft 
sub jet collinearly splitting at the jet boundary. Or more concisely, I consider edge of jet sub¬ 
jets. Having worked out the phase space region, I then give the appropriate factorization 
theorem for it. I then calculate to one-loop the objects found in the factorization theorem, 
finding both double logarithms of collinear origin, as well as a DGLAP style splitting 
process [54-56]. Having these tools in hand, I show how the single dressed gluon with these 
collinear effects can be included to resum NGLs of more inclusive observables, following 
the procedure of [42]. I then compare to the NLO BMS equation derived in [37]. I find that 
this evolution equation naively misses the fixed order subleading NGLs, which would 
be corrected with the appropriate NLO calculation for the trace of the density matrix. I 
find this to be connected with the lack of collinear evolution in the leading order BMS 
equation. Going back to the original BMS equation, and using the soft jet factorization 
theorem to calculate soft jet production at large-A^o I show how one can collinearly improve 
the equation to include these collinear double logarithms. Then I discuss whether one can 
find any other IR structure beyond out of jet fragmentation. Then I argue that the soft jet 
expansion gives a systematic view of calculating the resummed decay rates used in these 
universal IR evolution equations, highlighting its physically observable features. Finally, 
having started at the beginning, I then conclude at the end with speculations on small-x 
physics. 

2 Phase Space Region 

For the purpose of understanding the collinear effects in subleading NGLs, I consider 
e'’'e“ —)• hadrons with the hard scale Q^, and wish to isolate the region of phase space 
associated with a soft subjet with two hard jets. In particular, I suppose two hard jet 
regions, the “out-of-fat-jet” region where an inclusive jet shape measurement has been 
made, and a fat jet region containing (most) of the soft subjet. The two hard jet regions 
are back to back, with associated light cone vectors: 


Fat Jet Axis: n = (l,h). 

(2.1) 

Recoiling Jet Axis: h = (1, —h), 

(2.2) 

n • h = 2, = 0 . 

(2.3) 
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The soft subjet enjoys its own light-cone direction Ugj = {l,hsj) with conjugate Usj = 

(1, — figj). 

The phase space region for collinear splittings along the jet boundary exists strictly only 
for cone algorithms. A recombination algorithm will necessarily deform the jet boundary 
if a soft jet is present, though this effect can by mitigated with the anti-kj algorithm [57], 
where the factorization of [42] also applies. Here I suppose that the jet axis is defined using 
the broadening [58] or thrust axis [46] of the event, and a fixed cone of radius R is drawn 
about the jet (this is similar to the set up in Ref. [59]). For multi-jet events, an ideal cone 
algorithm for this factorization theorem is the XCone algorithm of Refs. [60, 61]. 

Within the fat jet, I measure the energy correlation functions and [47] 

to isolate the appropriate region of phase space^, whose dehnitions are summarized in 
App. B. In this cross-section, some out-of-fat-jet measurement has been made, r. I shall 
for dehniteness focus on the out-of-jet thrust measurement, with its action on a state jX) 
given by: 


(2-4) 

i&x 

The jet algorithm constraints on a single particle state are: 

©Fjjp) = 0fj{p)\p) = e(ia.v?^ - ^)1p) 

V 2 n■pJ 

^FJ\P) = ^FjiP) \P) = ^{- tan^ f + 1:^) 

The specific power counting I adopt to isolate the one soft jet region of phase space at 
the edge of the jet is: 


"3 


iA9 


SJ) 


4“^ ~ Qef'^ 

(2.6) 

F>y 

(2.7) 



3 1 

(2.8) 


AOgj = R — Ogj is the difference between the jet radius R and the angle of the soft jet 
axis hsj to the fat jet axis h, denoted Ogj- This necessitates the soft jet to collinearly split 
into and out of the fat jet, pictorially represented in Fig. 1. Specihcally, this implies the 
collinear modes of the soft jet in the soft jet factorization theorem of [42], and the boundary 

®For a complete and detailed discussion of the 1 —>■ 2 phase space for QCD radiation, and the impact 
the different regions have on the form of factorization, please see [62]. 
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Figure 1: A soft jet at the edge of a jet. The dotted lines correspond to the jet boundary, 
soft modes, cannot be factorized, since the two obey the same scalings: 


Psj ~ Ej 


Pbs ~ Ej- 
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( 2 . 10 ) 


where I have adopted the light cone coordinates of the soft subjet, {ngj ■ p,nsj ■ p,p±nsj), 
and Ej ~ Q is the energy scale of the fat jet. 

3 Factorization Theorem 

The factorization theorem takes the form: 
da 


de^^'^ de^^^dr 


aoH{Q^)H:^^ , ef , r) ; R 

® Qsj: ® (8) Jnir) ® Snn(e^^^;T; R 


(3.1) 

Where for conciseness, I have adopted the convention that any repeated argument is con¬ 
volved with all other functions sharing that argument (excepting the jet radius R): 

„ iV N 

Fi{x)® E2 {x)...® En{x) = Y[dxi6(x - '^Xi'jEi{xi)E2{x2)....FNixN) (3.2) 

i=l i=l 

A new feature of the edge of jet factorization theorem is the edge of jet function Susj- 
This is a novel jet function describing the fragmentation of partons into and out of the 
jet boundary. Though containing these fragmentation effects, the function is IR hnite and 
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calculable in perturbation theory, given the jet boundary cuts off the collinear divergence 
of the fragmentation. 

I have also explicitly indicated the dependence on the soft subjet large momentum 
fraction Qgj. Implicitly, this is a function of Q, e^\ and The proof of this factoriza¬ 
tion theorem, like that for PDFs or fragmentation functions [63], takes place at the level 
of the cross-section or squared amplitude, not the amplitude itself, as in previous subjet 
factorizations. This is since the hard matching for has both real and virtual correc¬ 
tions contributing to it. One constructs a complete basis of gauge invariant IR functions 
corresponding to all momentum regions. Since each sector has a unique leading power func¬ 
tion according to the power counting, after including appropriate subtractions to remove 
overlaps, one is led to (3.1). 


3.1 Definitions of Factorized Fnnctions 

I now give the operator dehnitions of the factorized low scale functions appearing in 
Eq. (3.1). Each measurement function should be expanded in the power counting of the 
factorization theorem. The expanded versions can be found in Ref. [62]. 


Tri-pole soft function: 

/ ic 

-i 


Snn 


4“^^) = 


(iepe®3 eg -- 


msn^.SnSnm 


— 200+C 


Q - f tA Q - ( 

OnnlCg , CL I I Cg , 


R 


5nn(4“^;^) = tr(0lr{5„5^}e 


Sn^,nJ~eP-.R] = 


-4“^0FjE3(“) 


sf{SnSnm 


it (<^) 

— 63 > 




(3.3) 

(3.4) 

(3.5) 


To avoid convolutions between the necessary subtractions to the naive function, I 
give the definitions in laplace space. I have explicitly given the appropriate soft 
subtractions that must be performed on the tri-pole soft function. One removes both 
the boundary soft radiation associated with the soft subjet, and the soft radiation 
from the original n-n dipole. 

• Edge-of-Jet Function: 

= (3.6) 

(0)5(r - Q^)5{Qsj - Qfju.j ■ V)5^^HVFsA5(et^ - ^c„,^./^(0)l0) 

The jet region constraints Qpj and Qyj restrictions project the state according to 
whether the momenta is inside the resolved fat jet, or the inclusive out-of-jet region, 
respectively, according to Eq. (2.5). The field operators are gauge invariant collinear 
gauge field strengths, whose definition can be found in Ref. [64]. 


- 5 - 





Figure 2: An example soft diagram that contributes to the soft jet production matching. 


• Dipole soft function: 

Snn{4°^;r-,R^ = tr(0|r{5„5^i}(5(e(“) - 0irjE3(“)|^)<5(r - 0j^T)r{5„5^}|O) 

(3.7) 


Note that this function generates soft contributions to both r and 63 °^^ 

The jet functions Jn and Jn are standard for the given jet shape, see for instance Refs. [12, 
13, 58]. 


4 Matching for Soft Jet Production 


Here I detail the one-loop matching procedure for the edge-of-jet factorization theorem^. 
The soft emissions within the fat jet are strongly ordered with respect to the soft subjet, 
thus the calculations are identical to those of [42] . The virtual corrections can be extracted 
from [67]. However, the “in-out” configuration exhibits new terms, since r ~ ~ ■ 

Therefore, the double real emission diagrams of full QCD are not expanded in a strongly 
ordered limit, but only in the soft limit compared to the hard scale Q^. The integrand 
Snn{p,Q) for the double soft real emission is given in [ 68 ]. This corresponds to two soft 
gluon exchanges between the wilson lines in the n and h directions, with Fig. 2 giving 
an example diagram. I can write the two real emission contribution to the fixed order 
cross-section as: 


(4.1) 




(4.2) 


Where and are the phase space constraints for when both emissions are inside 

the fat jet, and when one is inside the fat jet and the other out. “S.O.” denotes strong 

^Throughout these calculations, I have made use of the Mathematica package HypExp [65, 66]. 


- 6 - 






ordering. One emission will always be identified with the soft subjet, and to this end 
it is more convenient to trade the energy-energy correlation functions for the light-cone 
momentum fractions in the n-n coordinate system: 


u = n-psj, v = n-psj, (4.3) 

= {u + v)(2—^Y ^ ^ (4.4) 

With an identified soft subjet, this is always possible up to power corrections, with a simple 
jacobian factor. 

The matching itself will be determined by the difference between the full theory matrix 
elements and the effective theory. That is, the one-loop matching to is fixed by 
demanding: 


- v)5{T)5{ef^) + v) ^ 

dudvdeX 'dr ^ 


+ + 0{ai)... 




(4.5) 


with the tree-level matching being: 


M - 

^nn 


sjiT). . ^ 

^ ’ ^ 7r(uu)i+^r(l - e) 


OisCp 


l-|-€ 


n • n 


2i+^7rr(l - e) \n-psjPsj ■ n 


(4.6) 


The sum denotes the sum over one-loop matrix elements given by the infra-red 

functions of the factorization theorem (3.1). The one-loop hard matching I split into three 
contributions, in-in, in-out, and in-virtual: 


Hnn^\u,v,T) = H, 


- "nn^\u,v)6{T) 


+ Hnn^\u,v-T) ^ ^ , (4.7) 


in-out 


in-virt 


One emission must always be in the fat jet, since the soft subjet has been observed. 


4.1 In-In 


Restricting to the in-in configuration, in the strongly ordered expansion, the two real 
emission contribution is exactly reproduced by the factorization theorem contributions, so: 


M - 

^nn 


si(i) 


{u,v] 


= 0 


(4.8) 


This reflects the fact that the matching can be performed at the level of the amplitude, 
using only the IR finite terms of soft currents [42]. 


4.2 In-Virtual 

The real-virtual contribution is non-trivial. The infra-red divergent contributions are can¬ 
celed by the virtual loops of the effective theory matrix elements, so that this contribution 
is completely determined by the finite parts of the soft-virtual contributions of the full 
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theory. These are equivalent to the one-loop soft current, giving: 


Hnn^\u,v-n) 




[uvj 


\-i 


- vln^ 


in-virt 


TT^ 



(4.9) 


4.3 In-Out 

Here the full theory contribution is not fully reproduced by the effective theory, so that 
the difference between the effective theory and full gives a non-trivial contribution. I 
split the integrand into abelian and non-abelian components. The abelian contribution 
will be reproduced by the tree-level soft jet matching (4.6) and the out-of-fat-jet one-loop 
contribution from the dipole soft function (3.7): 


^^in-out, ab. 



(4.10) 


The non-abelian piece has a more subtle matching contribution. I write: 


^^in-out, n.ab. 


sj{T) 


M - 

^nn 


{u,v 


u + v;t] R) 


in-out 




in-out 


(4.11) 


Where I note that the soft jet large momentum fraction is given by tt -|- u = 2Qsj. Next I 
explicitly calculate the full theory. This is simply the integral: 

dcjin-o^h = N J[d^p]+[d'^q]+e(^'^ - - n • g) 

^(tan^^ - - n-p)6{v - n-p)Snn{p,q) (4.12) 

\ 1 n•p/ 

N = G^CFCAipe'f^^TrY (4.13) 


We are using the MS renormalization scheme, and •je is Euler’s gamma number. The 
integral measure can be written in terms of the relative angle between the transverse 
momenta components of p and q, the other angle being trivial. Using the on-shell conditions 
to fix the magnitude of the transverse momentum integrals in terms of the light cone 
fractions, I then have for integral measure: 


/ 


— tarn — \6[t — n ■ q]d[ tan^ — ~ ~ 


u — n ■ p)6{v 


[dM+[d^q]+o{ 


n ■ q 
n ■ q 

oo 


2-4-h2e^-2 

r(l-e)r(i-e) 


n ■ p) 

(4.14) 

(4.15) 


I have labeled the out-of-jet emission’s light-cone momentum fractions as x = n ■ q, and 
y = n ■ q. I now perform a further change of variables: 


X —)■ MX 


y —>• vxy 


(4.16) 
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Then with the further variable change, x —)■ x{vy)~^, I can evaluate the x,y, and (j) integrals 
as a series in e. Further, it is expedient to divide the in-out contribution into two pieces, 
isolating the collinear splitting scale 1 — ^cot^^ that will be removed by the effective theory 
subtractions. The final result is, dropping all terms that vanish as e —?• 0: 


da 


in-out, n.ab. , in-out, n.ab. , , in-out, n.ab. 

= «^split + «^no-split 


j in-out, n.ab. _ jj'Sj(T)/ \^sCa ( 'U'^^ 

^nn ^ 


^ split 


TTT 


)']p(t.,r)ln(l-5coe|) 


V 


'T + u cot^ ^ 


I uvT{v-T)coe^ / _ ^ r \ 

(u-|-r)^(r-|-ucot^=§) ^ Ca / 2{v + t) \r-I-cot^5 / 


Ca ^ 2{v + t) Vr-|-ucot^ — 


(4.17) 


_l„(l _ +dn(l - - l)| +o(£) (4.18) 

, in-out, n.ab. t Dt M f 

^^no-spUt = ’(U, < P{v: r)ln{^ 


+ 0[e) 

(4.19) 


I have explicitly factored out the all-orders in e tree-level hard matching. Where P{v,t) is 
given by: 


P{v,t) = 


(u^ + VT + r^)^ 

{v + tY 


Ca {v + tY 


(4.20) 


Note that P{v,t) is finite as r —)• 0, and the no-split contribution has no soft divergence 
as T —)> 0. It is the gluon’s 1 — >■ 2 splitting function with soft support removed. One 
can explicitly check taking the laplace/fourier transforms of u, r, and integrating over the 
allowed v, these terms with 7? = | are sufficient to reproduce the hemisphere NGLs of 
(7.1), and both the split and no-split terms contribute. 

To make the divergent behavior of this contribution to the cross-section transparent, 
I take the laplace transform of out of jet variable r in Eq. (4.18), whose conjugate variable 
I denote by r. Then I can simply expand in e, without making use of plus distributions^. 


7-'in-out, n.ab. 



dr e ^'^da 


in-out, n.ab. 
split 




TT \ V 2 


p 




dr 

T 


e-^^P{v,T) 


(4.21) 


1 i f 

^ ^ 2 I tan 2 # - ^ 



+ Plvf) > -h 0(e) 


(4.22) 

(4.23) 


®I factor out the all-orders in e tree level matching. The higher orders in e do not contribute to the 
matching, even though they naively give 0(e°) terms when multiplied against the terms of the one-loop 
correction. However, since all e divergences cancel when I sum over sectors, and the higher order terms in 
epsilon from the tree-level matching multiple all sectors identically, these spurious finite terms cancel out 
in the total cross-section. 
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One can work out the exact analytic expression for P in terms of incomplete gamma 
functions and exponentials, but the form is relatively unenlightening, other than it is a 
function of the product vf only. To complete the matching calculation, I need the edge-of- 
jet jet function, which is taken up in the next section. 

5 Edge-of-Jet Function 

I start this section with a general discussion of the phase space for the edge of jet function. 
Then I focus on the calculation for the contribution relevant for the matching and the 
resummation of NGLs, the in-out conhguration. The in-in configuration is relegated to an 
appendix. 

5.1 Tree Level and One-Loop Phase Space for Edge-of-Jet Function 

The tree level Edge-of-Jet function is calculated to be: 

J?) = 5{e‘f^)6{T){2TTf J [d‘^p]+5^‘^\p±)5{Qsj - fisj ■ p) 

= • ( 5 - 1 ) 

The transverse momentum constraint localizes the tree-level result to align with the ob¬ 
served soft jet axis, so that no contribution is given to r at tree-level. Moving to one loop, 
I have: 

=9^C^(2») J [d‘'ki]4d‘‘k2]+(i’““''{e'-p,Qsj,R,ki,k2)HT) 

+ 4."—(T. Q.j.R. ti. fc2)«(4“’)) (5-2) 

where the splitting function for an off-shell gluon into two partons with momenta ki, /c 2 is 
given as: 

Pgifisj ■ ki,nsj ■k2) _ 2 ‘ ‘ ^2 + {fisj ■ ^ 2 )^)^ 

ki ■ k2 fisj ■ ki fisj ■ k2{fisj ■ ki + fisj ■k2Yki- k2 

TfUf 2nsj • kifisj • ^2 - (1 - e)(% • ki + fisj • /c2)^ . . 

Ca (1 - e){nsj • ki + fisj ■ k2)'^ki • /c2 
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The two phase space restrictions break down as: 


i?, fci, fca) = 6{Qsj - n,, ■ h - Usj ■ k2)S I - N 


(a.) ( 2/ci • /t2 


/ Z/Ci • fe 2 \ 1 
Q Q * ^l^sj ' ^2 

4n L • k2± 


+ feO 0 (A«., + + 


sj • fii[n • Tisj 


$“(r,Q,„i?,A:i,fc2) = S{Qsj - n., • fci)< 5 (r - • k2)6^^-^Hki± + k2±) 

4n_L • fc 2 _L 
k2{n ■ UsjY 


9{M, 


4n L • /ci L 


'•sj 


ki{n ■ UsjY 




Usj ■ k2(n- risjY 
(5.4) 


(5.5) 


I-SJ 


All transverse projections _L are with respect to the local soft jet light-cone coordinates 
defined by Ugj and fisj. I have introduced the angle to the jet boundary IS.9sj'. 


, ^ oR n ■ Usi 

A6sj = tan^-^- (5-6) 

2 n ■ Ugj 

The 0-function constraints come from expanding the fat jet constraint Eq. (2.5), in the 
coordinate frame of the soft subjet when it is pressed against the boundary: 


^ / nR n-k 

0 tan^-—- 

V 2 n ■ k 


— 0 (AOgj -\- 


4n_L • A: 2 ± 
fisj ■ k2{n- UsjY 


-h ... 


(5.7) 


To avoid expansions in plus-distributions, I will henceforth work with the laplace transform 
of e‘f\ 4"^. 


5.2 In-Out Contribution 

The bare contribution is given by: 


41] (^; u, v; R^j 


in-out 


asCA sec 7re _ v / '^Y p f \ 

43~^e Tr^/‘^T{k - e) (v + TYKr'^itaiYR -v) 


Again, I want to take the laplace transform to extract the divergences. To do so I hrst 
isolate the soft contribution, by writing: 


Pg[v,T) = 2 


(v + tY 


VT 


(-1+P(u,r)) 


(5.9) 


So: 


41] ('t; w, v] R^ 


(XsCa 


sec vre 


in-out 


2^ ^^e7r^/'^r(k - e) T\T^(tan'^§ --) V 




OisCA 


^2 
sec vre 




u 


2i-2^e 7rV2r(f - e) r Vr2(tan2f - 




-^) P(u,t). (5.10) 
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Now I take the laplace transform, and expand in e, to achieve: 


42 (^; u, v] 

42 (^; 


div 

in-out 

fin 

in-out 


Inf(vf) 


Ca<^s 


TT 


1 1 1 
H-In 


fif 


4e2 ' 2e"""Vtan2f - ^ 

2 V 


a. 


:CAi 




TT 

+ 4n2 


-P(nf)ln 





fif 


tan^f - ^ 




cJt ~ T 

/ —e-^^P{v,T)hi- 

Jo r V 

TfUf 

Ca Jo r {v + tP 


The finite terms in momentum space are: 


(5.11) 


(5.12) 

(5.13) 

(5.14) 


42 ('^i 


R 


fin asCA\P{v,T) 


in-out 


TT 


l 


In 




r(tan2| - (() 



+ C^ 




tan2f - ^ 




The in-out contribution to the edge-of-jet function has several noteworthy features. Aside 
from containing double logarithmic divergences, it also features a single logarithmic diver¬ 
gence corresponding to the non-soft terms of DGLAP splitting function. At higher orders 
one can expect departures from DGLAP evolution. The appearance of DGLAP is not 
suprising, since essentially at this order in perturbation theory, the measurements imposed 
give the probability for a jet close to the boundary to fragment a parton out of the fat jet 
boundary. The IR divergence associated with this fragmentation is cut off by the angular 
distance of the soft jet to the fat jet boundary, so though the function will have an evolution 
similar to a fragmentation function, it is IR-finite and calculable in perturbation theory. 

It is natural that the double logarithmic divergences, here being related to the soft 
limit of the splitting functions, are controlled by an anomalous dimension proportional to 
the cusp to all orders in perturbation theory [69-72]. 


5.3 Edge-of-Jet Function Matching 

The edge-of-jet function admits an OPE onto the boundary soft function and a standard 
jet function, of the form: 




CnjT,R)Sn,^nJeP,R)JnJe^'> 




(5.16) 
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The tree-level matching is: 



(5.17) 


Since the boundary softs and standard jet function calculations reproduce eqns. (E.2) to 
(E.8), the divergent and finite contributions in laplace space to the matching coefficient to 
one-loop accuracy is then given by Eq. (5.11) and Eq. (5.12): 



42 

42 


div 

in-out 

fin 

in-out 


(5.18) 

(5.19) 


The renomalization group equation satisfied by this renormalized matching coefficient is 
to one loop: 


vr I Vtan^E _ ^ y u/ 


/r4lnC'n,,.(f,i?) =-— ^ 

dfj, ^ ^1 '"tarn# — V V 


(5.20) 


The anomalous dimensions of the other functions in Eq. (5.16) are given in Ref. [42], and 
one can construct the full RG equation of Sn^j from those results. 


6 Edge-of-Jet Matching 

I now have all the ingredients to finish the matching in Eq. (4.11). I use (4.21), (5.11), and 
(5.12) to compute: 




v 4 + 4 lnf+ ^P{vf) \ ( 6 . 1 ) 


in-out 
fin 
in-out 


{u,v,f) = (u, v) 

— Inf 


TT |4e2 2e ^Vtan^l V 2 e 

asCA j 1, 4') 

vr 4 2'''4an2|V4 


u \ / ^ \ 

4tan^^7 V utan^^v 48 




- .P^°(v;f) - Inf{vf) 


+ no split > . 


( 6 . 2 ) 


- 13 - 

















I also give the finite contribution in momentum space in terms of plus distributions. This 
allows us to show the structure of the “no-split” terms more transparently: 


H, 






fin 


in-out 






utan^=| 


viav?^ 


7r2 

+ l 6 


-I- ^/icot^^ 




T + U COt^=§ ' 


2(u + r)^''Vr + ucot2f/ 

uur(u-r)cot^f / _ ^ Tfnf \ 
(u-b r)^(T-|-ucot^^) ^ Ca ) 


+ Tfnf 


V^T \ 

(v + t)^ j 


(6.3) 


The plus distributions are dehned in App. A. Importantly, this contribution to the matching 
coefficient is manifestly finite as u —)• utan^^. That is, all sensitivity to the IR scales 
associated with the collinear splitting at the edge of the jet have been removed. 


7 One Soft Jet Contribution to NGLs with Collinear Effects 


7.1 Fixed Order Non-Global Logarithms 

For reference, I quote the fixed order result at for hemisphere NGLs, see Refs. [73, 74]. 
Dividing out the global contribution to the cross-section, the hemisphere NGL contribution 
for an initial dipole in the fundamental representation is: 


NGL ( 


cr — 


a 


TT 


U 




= -4Cf - TT^CAln^- + In- Ca 


vr^ 


12 


llvr^ - 3 - ISCs 
36 


-bTijn/- 


6 — dvr^' 


36 


(7.1) 


One can check by explicit calculation that the result for A^ = 4 SYM is simply the leading 
terms in transcendentality. The leading NGL at for non-hemispherical jets were obtained 
for a variety of algorithms in Ref. [75], and both leading and subleading were obtained for 
cone algorithms, similar to the ones considered in this paper, in Ref. [59]. 

It is important to note that the Cs subleading log is directly tied to the collinear double 
logarithms of the edge of jet function found in (5.12). In particularly, it arises in the full 
theory calculation of the NGLs integrating over the double logarithm at order e in the 
second line of Eq. (4.18). This fact will play an important role in the collinear evolution 
of the BMS kernel later. 


7.2 Dressing Gluons 

To elucidate the collinear splitting contributions to the NGLs, I following the procedure 
outlined in [42]. I re-associate functions in the factorization theorem (3.1) to obtain anoma¬ 
lous dimensions that allows me to resum the NGLs for a single soft jet emitted off of the 
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n-n dipole. I define; 


Wnn [Psj ]T]R)= Urn {psj i T; i?) £n,, (4“^) ‘S'n,,n,, (4“^ i ‘S'nnn,, 

63 ^0 

Gfinrisj ir 5 — ^'^sj ’ '^) 



(7.2) 

(7.3) 


These functions differ from those used in the original soft jet factorization to resum NGLs. 
Because of the modified power counting, both W and G can depend on the out of jet scale, 
since formally I have assumed no hierarchy. The RG equations satisfied by these functions 
are given by: 


In 


fire 


HE 




MCa f 


- + P{vf) 
v/ 


d. ^ m as{n)CAj, ( pfei’^ [u\ ~ 

R) =-—1‘"( (,atf f - J) V 


(7.4) 

(7.5) 


I can see how the RG-improvement of the WG product resums the collinear splittings: 


iPsj ■) R-) f^^Gnnnsj -^5 /^) — {Psj R-; PjGnnnsj ? R-) (M? 1 ^') 

(7.6) 

It is convenient to split the resummation factor into a double and single logarithmic con¬ 
tributions as: 


/■M A,,' , 


= Exp 


'Mi 


dM as{p')CA ^^/' p \ 
p' vr V//j/ 


— In 


U\n{P,Puf) = Exp 


[u\ dfi'as{p')CA\ 

— P{vf) 


(tan^l - ^) V M tt 7 

d^x'as{p')CA 




1^ 


TT 


(7.7) 


(7.8) 

(7.9) 


I can then write the cumulant momentum space resummed distribution as a convolution 
between the single and double logarithmic functions: 

roo 

U^{fi,^i;T)= dT'U^^2ip,Pi]P (7.10) 
Jo 
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I cannot obtain Uin{fi, Hi] t) in analytic form, but I can solve for U^^ 2 {h-: 




T{1 - ui{H,Hi)) 


r dn' 
= / — 


dp! as{H')CA 


TT 


QsC'A ln 

TT Hi 


V/r(tan2f \ vJ J 


u\ dn' Ois{H')CA 


TT V U,- / V 


//jC 


.7b 


Ai' 


TT 


VT V/r(tan2|-^) 


(7.11) 


(7.12) 


(7.13) 


The single soft jet contribution to the cumulative non-global logarithms then takes in 
resummed form: 


dv dT<Wnn{u,V]T] h) ® G{u,V]T] Hi) 

J U — T^ ^ 

® U^^ 2 {H,Hu'r) ® Uir,{H,HuT) - c-bin| (7.14) 

The limits of integration are determine as follows. The variable v is just integrated over 
all allowed angles of the fat jet, while u,t are bounded as follows: 



T > T 


T + U> 


U^>U (7.15) 

u> (7-16) 


The conditions of the first line are just the definition of the cumulative values of the 
measurement. The second line is more important. It states that the parent gluon that 
splits in and out of the jet is above the IR scale of the problem. Thus the edge of jet 
resummation can be applied without interfering with the global and IR divergences. Below 
the IR scale, and the gluon contributes to the global divergence, and no resummation 
should be applied. This also guarantees that all IR divergences cancel between the in — in, 
in — out, and out — out regions for the hemisphere soft function [73]. 

To see whether the dressed gluon will reproduce the NGLs, I must expand each of 
these functions to first nontrivial order, and evaluate the integrals. Note that the collinear 
bin is the expansion in the limit that u <C u, r before the integrals are evaluated [76]. This 
removes any overlap with a collinear subjet factorization. To this order (mixed leading- 
logarithmic prime, since we are incorporating matching contributions, and q;<j running, but 
ignoring the two-loop cusp contribution), the collinear-bin insures that the limit of Eq. (7.2) 
is well defined. Any part of the one-loop matrix elements contributing to the resolution 
variable will cancel against the same contribution with the collinear-bin expansion 
performed, up to corrections that are beyond the logarithmic order I am working. 
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7.2.1 Matrix Element Contributions to NGLs 

Due to the collinear bin, the only contribution to the NGLs that can arise from the matrix 
elements of the effective theory are those of either the in-out contribution to the soft jet 
production coefficient in Eq. (6.2), or the edge-of-jet matching coefficient of Eq. (5.19). 
Thus we have to order ag in each of the low scale functions: 


in-out 


Wnniu,v;T; fi)Gnnnsj{T] R; l^i) = v)6{t) + V, T] n) 

+ (r; R; m) + 0(af) 

Examining the matrix element in Eq. (5.15), one can choose the low scale to be: 


... 2 ^ u\ pv 

2 - JVn 


(7.17) 


(7.18) 


This will minimize all logs in the edge-of-jet function that can contribute to the NGLs. 
After the collinear bin, no terms from the edge-of-jet function will contribute to the 
NLGs. Examining the soft jet production matching coefficient, taking: 


/r ~ r 


V 9 i? 

tan — 
u 2 


(7.19) 


will minimize the logarithms in Eq. (6.3). However, with this choice, there are still terms 
left over that can contribute to the NGLs even after the collinear bin: 


Hpn^\u,v,T-,n) 


— c-bin = 


a: 


CfCa ( 


in-out 


TT^ UV 




-|- -P{v, T)ln 
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T + U COt^ Y 


— (5(T)ln( 
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V utan^ P 


In 1 - 


u 
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r -|- ucot^^ 


2(v + t) 

uur(u-r)cot^f / 
i?;-I--I-7/.cnt2|?i V Ca / 


(u-k r)3(r-|-ucot^^) V Ca 

In the hemisphere case, setting R = 7 r/ 2 , I evaluate the resulting integrals to get: 

'227r2-k9-k6 1n2-kl32 1n22 


(M (h r 

u V .L,_^C T 


J U—T 

n 


U + T/ V / 
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(hi ch 

r- u Ju V 
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72 

-In 
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In- + ... 

jC 
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r du dv r 

Jt- U Ju V 
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dr V — T 


f—= -^In— + ... 
\u + tJ 2 

ujTr 


_T-c T 2{u + t){t + vP 
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Ca 


= -1(i-2vTr 


(7.20) 


(7.21) 

(7.22) 
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1 u 

In-k ... 


(7.23) 
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Where we have dropped terms that are finite as ^ cx). Note that the large logarithms 
can be found simply by setting: 

ru 

/ dr —>■ — / dr + ... (7.24) 

J U—T^ J 0 

Then the finite terms are automatically dropped. This naturally coincides with an ordering 
of emissions u > t found in the evolution equation approach. 

7.2.2 Evolution Factors Contribution to NGLs 

Next I must expand the cumulant resummed distribution with scale choices given by 
Eqs. (7.18) and (7.19). I also take the collinear-bin and subtract. I find the following 
integrals to evaluate: 


Expanding (/i, /i*; Tc) - (it, u-Tc) - c-bin : 



(7.25) 


Expanding 17in(/r, r) — c-bin : 

r^du r dr f f N _ / 

J^C U V J^_^C T \ v) V ’ / ( 


llvr^ + 3 - 3 ln2 - 66 ln^2 


72 


rifTji + 6 — 6 ln2 — 24 ln^2 
72 


, u 

In-h ... 


(7.26) 


Again we have dropped terms that are finite as ^ oo. Upon summing Eqs. (7.21) 
through (7.26), and restoring the color factors and coupling, we reproduce Eq. (7.1). 


7.3 Alternative Resummation Schemes 

One can attempt to exponentiate more terms than is done with the canonical choice of 
factorization scale Eq. (7.19). In particular, one is tempted to set: 


^ ~ + n 


(7.27) 


This choice formally exponentiates all terms that can be connected to the splitting function, 
moving (7.21) into the single log resummation factor. However, the terms connected with 
the rational terms (7.23) would still not be exponentiated. Given that the logarithm 
associated with these terms is effectively a ratio of soft energy scales, r to r + rt, a more 
systematic route to the resummation of these terms would be to calculate and resum the 
subleading power corrections to the soft jet factorization of [42], along the lines of [64, 77- 
80]. 
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8 NLO BMS 


The comparison to the recently derived NLO BMS equation found in [37] is straightfor¬ 
ward®. Within this approach, one seeks to define an equation of motion for the reduced 
density matrix that results from integrating out the hard emissions. One does this time 
step by time step, where each previous step the mode that resulted from a decay into 
the IR now becomes a hard eikonal current, freely radiating. This results in a Markovian 
process at leading lo, and departures from leading log can be understood as corrections to 
the Markovian picture. The “time” evolution of the reduced density matrix is governed 
by a hamiltonian, and one can simply examine the action of the two loop hamiltonian 
generating the evolution on the initial dipole, perform the requisite IR averaging, and drop 
any terms that occur at order greater than a^. 

For simplicity, we restrict to N = 4 SYM, where the hamiltonian is simplest. If 


47r Ibvr^ 


( 8 . 1 ) 


is the hamiltonian generating the soft evolution, I have for this hamiltonian acting on an 
initial dipole in = 4 SYM: 

=2 [ "^4^° (217g°'f/or' - 


X (Tr 


[r“, T^] UiT^'t^' ui\ + Tr [r“', t^'pI 


+ 


Att'^Ca f ai2 


47r aioao2 


(tv[T^'UiT^uI]US°-' - CaTt[UiUI]^ ( 8 . 2 ) 


Where I have: 


^ab — 


13a ■ /3b 


/3a = (1,6a), l = 6a 


_(2)iV=4,/ 

^^nn-,00' 


Oi-n.n. 


« 0 na 00 '« 0 'fi 


2\n( 




\ Q^On^O'n 


+ 


1 + 




^On^O'n ^^O'n^On-' 


(8.3) 

(8.4) 


That is, all integrals are over null rays on the celestial sphere. I have gone ahead and set 
the initial dipole to be the back to back jets. The first term corresponds to multiple soft 
emissions generating new wilson lines Uo,Uo', while the second term corresponds to the 
cusp contribution. The IR averaging is accomplished by: 

t/“‘^r‘(^«(|)+e(cosfl-|)j (8,6) 

I focus on the multiple soft emission term, since this will naturally produce the single loga- 

®I thank Simon Caron-Huot for correspondence about the NLO procedure for BMS, as well as providing 
notes on the calculations. 
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rithmic contributions. Noting that one can drop the small cone about the other hemisphere, 
the integrations become with the IR averaging constraints: 



Here, x,y are related to the polar angles of the soft emissions. The NGL distribution is 
generated by acting Exp[—L atgl^R] on the initial hard configuration, so expanding the 
exponential and using the results in [37] for the kernel, we get: 







oICaCf 



( 8 . 8 ) 


When calculating with a hamiltonian germane to QCD, one obtains almost Eq. (7.1), 
except for the Cs term above, which has the wrong sign. In the full BMS equation, there is 
also the contribution coming from iteration of the leading order hamiltonian, however, this 
simply generates the leading log NGL at two loops. Then to get the <^3 logarithm correct, 
one must supplement the LO IR averaging procedure with an NLO averaging correction. 
However, we will see that collinearly improving the BMS equation to contain the leading 
double collinear logs results in a contribution that when added to the above, gives the 
correct (( 3 , allowing us to forgo the NLO averaging procedure^. This is not unexpected: 
under the assumption that the LO and NLO BMS kernels fully capture all soft coherence 
in the IR evolution, the only remaining corrections in the IR that can be obtained are 
incoherent emissions off of the wilson lines. 


9 Collinear Improvement of the BMS Eqnation 

It is beyond the scope of this paper to firmly establish how one might collinearly-improve 
the NLO BMS equation, however, noting the correspondence between the soft jet expansion 
in [42] and the leading order BMS equation, one can derive an expression that will hopefully 
capture all double collinear logs, at least in the large A^c limit®®. The purely non-global 

^One must be careful here since there is a great deal of scheme dependence as to where one puts the 
corrections for the NLLs, between the evolution hamiltonians, versus the IR calculation of the 

observable after evolution. Only the sum of the three are scheme invariant, and expected to give the full 
answer. 

®Since the collinear evolution is derived from color singlet objets, jet functions, this implies the collinear 
evolution does not care about the planar limit. However, the exact matching scale can be sensitive to the 
directions of the soft wilson lines the soft jet is entangled with. 

®The above NLO equation is for full color evolution, however, for checking the two loop NGLs, this is 
immaterial. 
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version of BMS is given as: 


, 9ab — 
dfi 

WabU) = 


OCsCa 




IT Jpj AtT 
a ■ b 


Wab 


9aj9jb 9ab 


}. 


a ■ rij rij • b 


(9.1) 

(9.2) 


Where a, b, and nj are null vectors. The integration is over the angular region of the fat jet. 
The RG equations of the soft jet factorization could be used to generate the Uabj factor: 


d , „ agCA 

fl — lnUabj = - 

vr 


j_'^{Waj{^ + Wj,{i)-Wab{i)} =-^^Sab{j). (9.3) 


Where the integration is over the angular region of out of jet region. From Ref. [23], the 
function Sab{j) always contains the term In^l — . Indeed, examining Eqs. (7.4) and (7.5) 

and noting that for the case of a = n, 6 = h: 

- Wnnii)] = -ln(l - 
u n ■ rij 

V n ■ rij 

I will show how to collinearly improve the BMS equation, assuming that the effective 
theories correctly calculate the resummed production rate of soft jets modulo global effects. 
This resummation of the collinear double logs maintains the conformal invariance of the 
BMS equation, up to the boundary condition of the jet radius. I construct the resummed 
weight in the BMS equation through a sequence of effective theories under the assumption: 

7/ 

r"(l--)<r", (9.6) 

V 

(9.7) 


(9.4) 

(9.5) 


I evolve between the scales of Eq. (9.6) using the evolution the logarithmic terms of 
Eq. (7.5), then I evolve between the scales of Eq. (9.7) using Eqs. (9.4) and (9.4). This se¬ 
quence of evolutions can be codified into a limit of the factorization theorem in [42] , where 
one essentially takes the threshold limit of the collinear splittings at the jet boundary, 
which is given in App. D. Then we have for the resummation in the BMS equation: 


Unkj = Exp 
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Then the collinearly improved BMS equation for the n — n dipole is: 
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(9.8) 


(9.9) 


- 21 













Finally, one would like to collinearly improve the BMS equation for an arbitrary dipole 
with directions a,b, derived in App. D: 


^abj 


= Exp 


OCsCa 


TT 


■Sa6(i)ln( ^ 


+ —Sotondi-- 


vr 


+ 


CX-sCa 

27r 


SlU) (9.10) 


One can check from the explicit expressions for Sab{j) (see App. C) this scheme also expo¬ 
nentiates the same double collinear log. Expanding out this scheme to find the predicted 
logs, one finds it over-estimates the Cs NGL by a factor of two. Namely, expanding out the 
collinear improved BMS equation and integrating, we find the NGL contribution: 


o^sCa du dv/ asCA -^/^ 

27r J^C u V \ TT V V/ 

^ oxICa ( 

27r2 V 

Correcting for the large Nc limit by swapping C\ —)■ 2CfCa, and adding this to Eq. (8.8), 
I find that the leading in transcendantality terms of Eq. (7.1) (the A” = 4 SYM result) are 
reproduced. The specihc to QCD contributions are already captured in the formalism of 
[37]. Thus the collinear improvement of the BMS equation removes the need to supple¬ 
ment the NLL resummation of NGLs with IR averaging corrections, that is, no large log 
still resides in the IR after evolution has taken place. The scheme dependence of the IR 
averaging has been completely shuffled into a redefinition of the BMS kernels, though there 
will now exist a scheme dependence between the collinear evolution, and the fixed order 
expression for the BMS kernels. The collinear resummation scales chosen here (the natu¬ 
ral ones from the point of view of the factorization theorem) correspond to the “Lorentz” 
scheme evolution kernel 

Here we have only used the leading singular pieces of the splitting function for the 
collinear improvement, however, other than having to solve a more complicated evolution 
equation, the full edge of jet function contains all the ingredients necessary to also capture 
these effects. Careful attention should be paid to the resulting changes in the kernel 
if further subleading collinear resummation is performed. 


2Tr \ V 

12 A ^ 


(9.11) 


10 Collinear Matching for Deep versns Edge Snbjets, and Finding the 
Buffer Region 

I seek to argue that the collinear corrections to the BMS equation for jets deep in the fat 
jet are irrelavant, exactly so in a conformal theory like N = 4 SYM. Every wilson line 
generate by the BMS evolution naturally comes with its own jet function, if for no other 
reason than to absorb and cancel collinear divergences. The question then is whether these 
jet functions at all have a natural IR scale beyond that of the fat jet’s mass (which sets 
their energy scale), or are the “inclusive-unobserved” jet functions of Ref. [13]. Given no 
natural IR scale, any loop integral in these jet functions will be scaleless, and in dimension 
regularization, they will be zero. In the calculation of the BMS kernel, they will appear 
as at most collinear corrections that cancel divergences in “hard” virtual corrections of 
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these kernels. In QCD, there will always be the IR scale of Aqqd, so though the collinear 
integrals will be scaleless in perturbation theory, running to non-perturbative scales is 
still necessary. Angles of subjets with respect to each other are always bounded by the 
jet radius, with no enhancing singularity, thus never give a meaningful scale. As subjets 
approach each other, they are simply combined, and collinear divergences cancel with no 
left over large logarithm. Then only close to the fat jet boundary does a new IR scale arises: 
the angle to this boundary that cuts off fragmentation into or out of the jet. This now 
gives the collinear corrections a non-trivial structure, even in a conformal theory, and one 
can conjecture that all scheme dependence of the IR averaging will arise from corrections 
having to do with edge of jet physics, since that is the only large scale left in the problem. 
Thus with the appropriate collinear resummation of the BMS kernel, one could render 
the IR averaging procedure trivial, in the sense that it is always given by the leading 
order procedure, or the leading order procedure times a kinematically trivial series in a*. 
Moreover, these corrections are universally predicted from a single jet function calculation. 
From an effective theory viewpoint, this is pleasing, since one would like to obtain all large 
logs from evolution, with no surprises in the IR trace of the reduced density matrix. 

Finally, it is important to quantify how big the edge of the jet is. Other than the 
potential confinement scale, only the size of the NGL can matter. So we have the different 
angular regimes: 

u u . 

I--< — orl-->— lO.I 

V fj, V IX 

In the first, the second term of Eq. (9.8) dominates the exponent, while in the latter the 
standard term dominates. Thus the size of the non-global log effectively sets the size of the 
boundary layer containing the collinear theory. The larger the NGL, the less the collinear 
resummation contributes to the improved equation. This is just the statement that the 
leading logs are genuinely leading. However, at small to moderate values of the NGL, the 
collinear effects could be expected to be substantial. One can get a feel for this competition 
by plotting out the size of the buffer region. This is defined to be the rapidity of the soft 
jet at which the resummation factor U in Eq. (9.1) attains its half maximum, given the 
size of the NGL. This is plotted in Fig. 3. 

11 The Lindblad Equation and the Expansion in Soft Jets 

These non-linear IR/RG evolution equations found in jet/small-x physics all fit a very 
general pattern (see also [81-83]). Namely, since the underlying physics is that of a Markov 
process at leading log, the equations are examples of the Lindblad equation Refs. [84-86]. 
These evolution equations for the reduced density matrix in the presence of an energetic 
enviroment decaying into the IR take a universal form, and in a field theory context have 
recently been discussed in Refs. [87-89]. If one takes a generic unitary quantum system 
with some initial energetic state with density matrix p, and begins to integrate out the 
energetic modes of the matrix, giving a reduced density matrix pn, then evolution in the 
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Figure 3: The rapidity of the soft jet at which the resummation factor in Eq. (9.1) reaches 
half its maximum value. Plotted are both the original BMS resummation, and the collinear 
improvement for the n-h dipole. Note: ~ .11. 


“time” t (or whatever observable the reduced density matrix is to be ordered in) satisfies: 

dtPR = [iH,pR] + ^^Tij{t)^2LipRVj — V-tipR — p^V-Li^ (H-l) 

ij 

if and only if the system is Markovian, and the evolution satisfies a semi-group composition 
rule while preserving the positivity of the reduced density matrix. H is a hermitian opera¬ 
tor, and the evolution is non-unitary (the flow is diffusive and potentially irreversible, thus 
the semi-group composition rule of the evolution.), and this is in contrast to traditional RG 
equations found in factorization theorems where one can flow in either direction through 
RG consistency. The Tij are the decay rates of the energetic enviroment into the IR, and 
Li create/annihilate the asymptotic IR states. The sum over i,j is over the kinematically 
allowed phase space for the decay, assuming a mode can be created at each point. In LO 
BMS, the creation and annihilation operators map to wilson lines, while the eikonal factor 
corresponds to the decay rates from high energy modes to softer modes. The difficult part 
though is not finding the form of the equation, but demonstrating factorization, that is, 
the system is Markovian to leading log. Effective field theory and in particular the soft 
jet expansion of Ref. [42] give a systematic way to calculate the resummed decay rates of 
the Lindblad equation for jet physics (the BMS equation) and its subleading logarithmic 
extensions such that the sum over the allowed phase space for the decays is well-behaved 
in all corners. 

For jet physics, the leading log factorization was accomplished long ago with the “jet 
calculus” of [90-92], however, instead of writing evolution equations for the reduced density 
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matrix^*^, the authors sparked the now well-established tradition of Monte-Carlo simulta- 
tion of jets that properly accounts for soft coherence. These Monte-Carlos solve these 
evolution equations [81], at least approximately, ultimately the more important endeavor. 

12 Conclusions 

I have presented a factorization theorem that can be exploited to resum terms contributing 
to subleading NGLs in jet cross-sections. Moreover the factorization theorem itself will be 
helpful in understanding the boundary structure of cone algorithms, and have useful appli¬ 
cations in jet substrucutre studies. The resummation is distinctly collinear, and captures 
the fragmentation effects of soft jets at the jet boundary. Using the factorization theorem 
to calculate soft jet production, I gave a simple extension of the BMS equation that would 
include the leading collinear effects, but is still not single log accurate. Interestingly, it 
captures logarithmic contributions that are naively missed by the NLO BMS kernel, but 
would be corrected when calculating the IR trace of the reduced density matrix as gener¬ 
ated by BMS evolution at NLO order. Taken together, the NLO kernel and the collinear 
resummation, they form a complete NLL exponentiation of the subleading NGLs, with 
a trivial structure in the IR matrix element. The soft jet expansion then, which can be 
systematic up to the number of jets included, provides an important check on what effects 
are captured in these evolution equations, and can be used to resum the evolution kernels 
thesmelves. 

With regards to phenomenology, it is plausible that the bulk of these collinear effects at 
large-Nc are already captured by antennae-dipole showers, like ViNCiA [94, 95] or Ariadne 
[96]. Given these showers capture leading log soft coherence and parton splitting at each 
splitting/emission, for all phase space points, they would naturally incorporate many of the 
subleading logarithms, especially the double collinear logs. However, some terms present in 
QGD, like the subleading NGLs with rational coefficients, would only be included once the 
parton shower was matched to a fixed order calculation with sufficient number of loops, or 
if the shower was constructed around a 2 —>• 4 splitting scheme. Indeed, this suggests that 
the appropriate construction of NLO BMS equation should be around antennae-multipoles 
instead of soft currents [97, 98] , to control the decay rates into IR physics. These naturally 
include both soft coherence and collinear splittings. Then one would be constructing the 
shower equation [81, 84-86] for a full scale NLL parton shower. 

Perhaps more important is the implications for the B-JIMWLK hierarchy [99-110]^^. 
NGLs and small-x physics are known to be related via a conformal mapping [32, 33], 
including the NLO corrections to B-JIMWLK, [37], whose conformal properties are more 
delicate [115]. This relationship is so-far exact in A = 4 SYM, and departures from it 
are directly related to the beta function of QGD. However, both the BFKL [116, 117] 
and the B-JIMWLK hierarchy resummations are noticeably improved at NLO with the 
inclusion of at least a partial collinear resummation of the evolution kernels [118-120]. 

^°But there was a consideration for time-evolution of the wavefunction [93]. 

^^For a detailed discussion of the assumptions behind B-JIMWLK, see [111], and for the NLO corrections, 
see [112-114]. 
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It would be pleasing to exhibit a conformal mapping between the collinear resummation 
applied here, and the small-x regime, at least for N = A SYM. For the NGL case, the 
jet radius provided a natural collinear cut-off for these subleading logarithms. Under the 
stereographic projection that relates the jet physics to small-x, the jet boundary is mapped 
to a circle in the impact parameter plane. Thus the concern would be collinear splittings 
between large and small transverse momentum scales (small and large values of the impact 
parameter). However, in the B-JIMWLK case there is no emergent IR cutoff like that from 
the jet boundary. In a conformal theory, all subjets produces by B-JIMWLK are then 
deep in the “fat jet”, and the collinear improvement seems pointless since there is no IR 
scale to run to. In QCD, one must worry about Aqcd giving the IR cutoff to the impact 
parameter space. Given the importance of these collinear improvements, this suggests 
complete quantification of all participating momentum regions in the high-energy forward 
scattering regime is lacking. 

Finally, given recent work on resumming all leading logs at finite-W in the hemisphere 
case [35], the stage seems set to numerically assess the impact of the NLL NGLs, at least 
out to moderate values of the non-global logs ^^^In)))^ ~ 2. It would be relatively 
straightforward to work to NLL with two soft subjet contributions, including collinear 
resummations, and comparing against the collinearly improved LOj-NLO BMS, which 
seems to be adaptable to the methods of [35]. Like in the small-x case where the collinear 
improvement allowed the NLO BK evolution to maintain positivy. Ref. [119, 120], it would 
be interesting to see if the same obtains in the jet case. 
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A Plus Distributions 


I briefly summarize the properties of the plus distributions found in the renormalized 
momentum space functions. I follow the definition given in [121]. Let q{x) be a function 
less singular than Xjx^ at the origin, then the rgularized plus distribution with boundary 
point xo is given by: 


q{x) 


[a:o] d 

= hm — 
+ s-^o dx 


9{x 


px 

Q(x,xo)= / dxq{x') 
J XQ 


5)Q{x,xq) 


(A.l) 

(A.2) 


- 26 - 







Integrating against a function f{x), the above plus distribution yields: 


c'f(x') q{x') 


, J [^o] 


dx q{x )[f{x ) - /(O)] + fiO)Q{Xmax,Xo) (A.3) 


If xo = 1, then I adopt the convention: 


L J + L J - 


Now I can define the distributions: 


^ ^ Lxi-“J + 

jn 

£„(x) = ^C%x) 

^ ’ da^ ^ ^ a=0 Lx 

Ca{X,x) = 

Cn{X,x) = 

These last distributions satisfy the rescaling identity: 


-In’^x 

-X 


£“(A,x) = A-“£“(x) + (I(x) 


A-“ - 1 


CniX,x) = nCkln>^[xYCn-kix) + 6{x) Yi 


(A.IO) 


This allows us to take the logarithmic derivative: 

^^j^£n(A,x) = ^-A:„C'fcln^"^(^A“^^£„_fc(x) - (I(x)ln”(^A"^^ (A.ll) 


-£o(A,x) = -5(x) 

-£i(A, x) = (5(x)lnA — C{){\^ x) 


(A.12) 

(A.13) 


Finally, the plus distributions encountered in this paper have the following laplace trans¬ 
forms: 


rOQ \ —a 


dx£“(A,x)e ""= / dx^ e"^-l +/ dx^^ + 


Jo 3^"“" ^ 


^oo ;^-a _)^-a _ 


ll X^ 


J dxCo{X,x)e —In^rAe'*'®^ 

J dxCi{X, x)e~^'^ = -In^ ^rXe^^'j + 


(A.14) 

(A.15) 


(A.16) 
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B Energy Correlation Functions 

The energy correlation functions are defined as: 



j^2 ^ EiEj 

i<j&J 


/ 2pi • pj 

V EiEj 


a/2 


1 


EiEjEk 

i<j<k&J 


f 2pi • Pj 2pi • pk 2pj • pk 

\ EiEj EiEk EjEk 


a/2 


(B.l) 


Here J denotes the jet, pi and are the four momentum and energy of particle i in the jet. 


C Out of Jet Anomalous Dimensions For BMS 


I give the explicit forms in the hemisphere case of the out-of-jet anomalous dimension of 
Eq. (9.3), as calculated in [23]. If a, h are light-like vectors within the fat jet, we have: 


Sab{j) = —Infcos 9^ -In ( j ) 

V V 2 \2[aj]\jb]J 


r _ ‘) 9 j\ / 1 -t- COS 0 ,- 

= -hi( 1 - tan^^j - ln(^-^^ 


2 \2lajm 


When b = n, we have: 


5a.(2) = -ln(cos0,)-lln(^i^ 

' 1 -|- cos 6j 


= —Inf 1 — tan^^ ) — In 


2 \2lajm 


Note: tan^Y = ^ v- ^ define the angular products: 


(ab) = 1 — cosOab = 1 — cosOaCosOb — cos{4>a — i?^fe)sin0asin0a 
[ab] = (ab) = 1 -|- cos9aCos9b — cos{(j)a — </>b)sin0aSin0a 


(C.l) 

(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C.6) 


9ab is the angle between the null vectors, and a is the reflection of the spatial components of 
a into the other hemisphere. We have given the explicit forms in the spherical coordinates 
defined by the fat jet axis n. 


D Joint Resummation of Soft Jets 

To incorporate the collinear double logs into the BMS equation (Sec. 9), we need to combine 
the factorization of this paper, with that of Ref. [42]. Thus we want a power counting 
scheme: 

u 3> r S> rfl — —cot^—) . (D-1) 

V V 2 / 
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First we recall the factorization of Ref. [42]: 


da 


de^2^ de^2^ de^^'^dr 


= aoH{Q‘^)H:^^ ( , ef ^, r) Snnn,, (; r; R) ® (e'f’, R) 




4“) 


Jnsjiei'"\Qsj) (8) ® Jnir) 




(D.2) 


In the region of Eq. (D.l) of phase space, the soft jet factorization theorem takes the form: 


da 


(a) W = <^oH{Q^)Kl (ei^^ , ef , r j; r; Rj(r, R) 

^^2 ^^2 ^^3 da 


4 «). 


Sn^jnsjie^^\ R) ® £n,j{e^i‘\Qsj) ® -^71(63"^) 0 Jnir) 0 Rnn (^£3"^; r; R 

(D.3) 


The matching coefficient R) is related to that of Sec. 5.3 with the added expansion 

V ^ T. This up to constants drops all terms but the double logs in Eq. (5.12). We have 
converted the soft function into a hard function Snnusj by subtracting the contribution of 

Snnn.j (4“^; -R) = Hnnu^j (4“V; C'n^" (E R) ® Snn{ef^-,T\ R^ + ... (D.4) 

This localizes the scale 1 — -cot^^ into the function and the scale r into the hard 

matching coefficent. The new three wilson line soft function is equivalent that found in 
Eq. (3.3). To get the resummation of the double collinear logs, we evolve C^'\ between the 
scales: 

—-cot^^^ —)■ r'^. (D.5) 

To leading log in the double collinear logs, we can simply exponentiate the logs, but at 
next-to-leading log and beyond, a more complicated resummation factor like Eq. (7.11) 
may be necessary in general. However, in a conformal theory, performing the resumma¬ 
tion in laplace space with canonical scale setting, the dependence on the out-of-jet scale 
vanishes completely in the resummation, so that the exponentiated form of the double logs 
would hold to all orders, with only corrections to the anomalous dimension. To finish the 
resummation, we run Hnnusj together to the scale /r from r'^. From the matching 

equation (D.4), and appropriately reassociating the functions of the factorization theorem 
to cancel global effects, this is just the running of the Gnun^j factor of Ref. [42]. Together, 
these resummations give the form Eq. (9.8). 

Explicitly, the evolution proceeds as follows. First I take the matching for arbitrary 
dipole o, b radiating a soft jet: 

5a6n.,(ei“V;7?) =77a6n.,(e^“V;7?)c);^J(r,R)®5a6(ei“V;7?) +••• (D.6) 
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Take the edge of jet anomalous dimension at threshold: 


d 


//—lnC“;(f,i?)5„ 


as{fi)CA 


TT 


(ln(^WO - ln(l - 


(D.7) 


I evolve between the scales: 


SJ 


n ■ n 


SJ 


ni?\ ./, n • Usi O-R , 

cor- W 1-r^cor- |e 

2 / V n • Usj 2 


^abi'^sj) 


(D.8) 


Though awkward, from the explicit expressions for the out of jet anomalous dimensions, 
Sab{nsj) in Sec. C, this last scale does not actually depend on the angle to the jet bound¬ 
ary. This intermediate scale is chosen as the scale where the out-of-jet component of the 
anomalous dimension of Habnsj vanishes in (D.6), so the matching logs are minimized. 
I then evolve using the standard dressed gluon anomalous dimension (9.3) between the 
scales: 

> /i (D.9) 


I, up to running coupling effects that are easily included, obtain Eq. (9.10). The appearance 
of ” ^ cot^^^ gSabi^sj) intermediate scale is just a consequence of the function 

Habusj ™ (D.6) being sensitive to positions of the directions a, b within the fat jet. 


E In-In Contribution To Edge-of-Jet Function 

The in-in contribution can be written in terms of a standard jet function contribution with 
no jet boundary contributions, and the jet boundary contributions: 



) 'T'l Qsj : R 



(E.l) 
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The standard contribution is given by: 
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(E.4) 

(E.5) 


The boundary dependent pieces are given as: 




Ca^^s 


+ 


TT I 26^(1 — a) e(l — 
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S3 


These boundary dependent pieces have hnite corrections that vanish as —)• 0. The 

leading terms given here of the boundary dependent contribution can be factored into a 
boundary soft function. 
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